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Abstract. We create displaced number states, which are non-classical generalizations 
of coherent states, of a vibrational mode of a single trapped ion. The creation of these 
states is accomplished by a combination of optical and electrical manipulation of the 
ion. A number state is first prepared by laser-driven climbing of the Jaynes-Cummings 
ladder, followed by displacement created by a sudden shift of the electrostatic trapping 
potential. Number states n =0,1 and 2 are prepared, and displacement amplitudes of 
up to a « 2.8 are reached. The states are analyzed tomographically, and we find a 
good agreement with the theoretically expected results for displaced number states. 
Quantum features elucidating the concept of interference in phase space are clearly 
demonstrated experimentally. 
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1. Introduction 

In the last decades, it has become possible to perform experiments on well isolated 
simple quantum systems. For example, a trapped single ion can be manipulated with 
laser interactions such that non-classical states of the ion's internal electronic and 
external vibrational states can be created and probed. Such experiments allow for 
exploring the fascinating features of quantum entanglement [U |2], and have set the 
basis for emerging quantum technologies, e.g. scalable quantum processors using ions 
in a segmented ions trap device [3]. Engineering and probing non-classical states of 
quantum systems, besides being of interest in its own right, has many prospects and 
applications. Techniques for quantum control and measurement can be established 
and benchmarked, decoherence properties of different quantum states can be compared 
[I], and theoretical concepts such as negativity of quasi probability distributions can 
be studied. Furthermore, nonclassical states set the ground for quantum simulation 
experiments E] and offer the possibility of enhanced interferometric resolution [7], 
providing the basis for future quantum sensors. 

An example of a nonclassical state is the displaced number state (DNS), which is a 
generalization of the well-known Glauber coherent state. DNS can be represented as 
\a,n) in Dirac notation, where n is the quantum number of an energy eigenstate of a 
harmonic oscillator, which is additionally displaced by an amplitude a in phase space. 
DNS have been studied theoretically [HI El CHI fTTj . they provide a versatile tool e.g. 
for the calculation of quasi-probability distributions, as their coherent or number state 
nature can be exploited depending on the context. Furthermore, DNS provide a clear 
illustration of the concept of interference in phase space p2], e.g. yielding a semiclassical 
explanation of molecular spectra. Experimental work on DNS done so far with photons 
includes generation and characterization of the state \a, 1) with displacements of up to 
a ~ 2.4 [13J, which revealed clear deviations from the ordinary coherent state \a, 0), and 
signatures of interference in phase space have been observed. A more recent experiment 
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with photons showed the non-classical decrease of the overlap integral |(l|a, 1)| 2 with 
increasing a [H]. Further experimental work established an analogy between DNS and 
classical light propagation in Glauber- Fock waveguides [15j HH] • 

For trapped ions, many fascinating experiments have been performed on the generation 
and analysis of nonclassical states in the last two decades, prominent examples being the 
creation of Fock, coherent and squeezed states [17], the tomographic reconstruction of 
quasiprobability distributions [TS] and more recently the entanglement of the motional 
modes of two trapped ions [19J. Our work is extending the previous work on nonclassical 
states of trapped ions by adding DNS to the collection of states generated and analyzed 
so far. While photonic experiments have been limited to states \a, 1), we prepare states 
of up to n = 2, and we are able to see characteristic higher-order interference minima in 
phase space. In addition, we demonstrate that both optical and non-resonant electrical 
forces can be applied together with sufficient precision to allow for quantum state 
engineering of the motional degree of freedom. We cannot perform a direct phonon 
number measurement, compared to the direct photon number detection in Ref. [H]. 
Thus we rather reconstruct the generated state of the motional mode of a single trapped 
ion by a tomographic method, which features a high precision and allows for a rather 
precise matching of the obtained results to the theoretical prediction. 

2. Experimental scheme 

2.1. Preparation of displaced Fock states 




Figure 1. State preparation scheme: a) shows how controlled displacement is achieved 
by switching the voltage of the trap segments. The ion is initially kept at segment A, 
while a voltage pulse on segment B temporally shifts the potential minimum, b) 
illustrates the preparation of number states by alternating 7r-pulses on blue sideband 
and carrier transitions. 



To create DNS of the axial mode of vibration of a single trapped ion, we combine 
different methods for the quantum control of motional and internal degrees of freedom. 
We first create number states by a sequence of laser pulses consisting of 7r pulses on the 
blue motional sideband, each adding a phonon to the motional excitation, and 7r pulses 
on the carrier transition to reset the internal state. Electrical forces generated with 
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control voltages at the trap segments are then used to displace the harmonic oscillator 
by a well defined amplitude. The scheme is sketched in Fig. [T] The displacement 
operation is accomplished by the application of a voltage kick generated by a fast 
arbitrary waveform generator. Such hardware technology with sufficient precision and 
low noise level became available only recently, enabling the fast transport of trapped 
ion over long distances [201 EI] ■ We note that displacement can also be accomplished 
by different means. One possibility is the application of optical forces [Tj, here the 
magnitude of achievable displacement is limited by the requirement that the oscillation 
amplitude in position space has to be well below the wavelength of the optical driving 
field [22J. Another possibility is resonant electrical driving [TS] . 

The detailed experimental scheme is as follows: A single 40 Ca + ion is stored in a 
microstructured segmented Paul trap [23] and laser cooled on the 5*1/2 ^ P1/2 cycling 
transition near 397 nm. The mj = +1/2 and mj = — 1/2 sub levels of the S1/2 ground 
state are henceforth termed |t) and \\), respectively. These levels are split in energy 
by an external magnetic field, leading to a Zeeman splitting of about uz ~ 27T-18 MHz. 
Coherent transitions between these states are driven via stimulated Raman transitions 
[24J. These transitions are mediated via the P\ii state and driven by two optical 
fields near 397 nm, red detuned from the cycling transition by a frequency of about 
A# rj 27T-40 GHz. Coupling to the axial mode of vibration takes place with a Lamb- 
Dicke factor of about 77 ~ 0.21 [25J, allowing for the coherent driving of motional 
sidebands. 

For each measurement, the axial mode of vibration at a frequency of u ax ~ 27T-1.35 MHz 
is cooled close to the motional ground state with pulsed Raman sideband cooling on a 
stimulated Raman transition, reaching a mean phonon number as low as n < 0.1. The 
internal state is then initialized to |t) by optical pumping. We then prepare a Fock state 
\n) by performing a 7r-pulse on the blue motional sideband of the stimulated Raman 
transition, changing the state according to 

|t,m)-H+,™ + l>. (1) 

The internal state is then reset by a 7r-pulse on the carrier transition, 

|4,,m + l)-Ht,m+l>. (2) 

This sequence is repeated n times. For each Fock state preparation step m, the pulse 
areas for the sideband and carrier pulses are calibrated by monitoring the respective 
Rabi oscillations. The preparation is concluded by additional optical pumping to ||) 
with a circularly polarized beam at the cycling transition. The fidelity of the Fock 
state preparation is limited by the imperfect preparation of the motional ground state, 
leading to a fidelity deterioration for higher n , since the transfer pulse areas depend on 
m and are therefore not well defined. This results in a spread of the population over the 
state ladder. The concluding optical pumping facilitates the evaluation of the data by 
setting a well defined internal state of the ion, which compensates for dispersion-induced 
errors of pulse areas of the preparation pulses. However, deviations of the desired final 
motional quantum number are not compensated, this error is even slightly increased as 
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there is a finite probability of changing the motional state in an excitation and decay 
cycle of the pumping. 

The displacement is then performed by changing the voltage levels on the trap segment 
and its neighboring segment, further details are explained in Sec. |2.3| The voltages 
are switched back to the initial values after 400 ns. We then apply an analysis pulse 
of variable duration on the stimulated Raman transition and its respective blue and 



red sidebands. The analysis procedure is further detailed in Sec. 2.2 Finally, the 
population in \\) is selectively transferred to the metastable -D5/2 state by electron 
shelving with laser light near 729 nm. Upon subsequent saturated driving of the cycling 
transition, fluorescence will be observed if the ion is found to be in the \\) state before 
the shelving, and no fluorescence will be observed otherwise, such that a spin-selective 
readout is performed. The entire sequence is repeated 200 times to obtain an estimate 
of the probability P-j- to be in the \\) state. This measurement is carried out for various 
analysis pulse durations and analysis pulse detunings, corresponding to resonant driving 
of red sideband, carrier and blue sideband. 



2.2. State analysis 

We describe the motional state of the ion, which is not necessarily a pure state, by the 
density matrix p. We assume preparation of the internal state |t), thus the total state 
is given by \\) ® p prior to the analysis pulse. The final probability to find the ion in 
the |t) state is then given by the Rabi oscillation signal 

P%An{0) = ^Eft(l + COs(M fe , An #)) , (3) 
1 k=0 

where 6 indicates the analysis pulse area and pk = Tr (p \ k) (k\) is the phonon probability 
distribution (PPD). This is an incoherent superposition of Rabi oscillations starting from 
different number states k, where the Rabi frequency is altered by the matrix element 
Mk,A n [22] in a nonlinear way. This opens up the possibility to infer p/, from the measured 
signal. The reconstruction of the PPD from a set of measured functions, P\{0), with 
5n = -1,0,1, is done with a maximum likelihood parameter estimation, where the pk 
values for k < /c max = 6 are varied freely, besides the requirement of a normalized PPD. 
Further parameters such as the bare Rabi frequency and the readout fidelity are included 
in the reconstruction to suppress systematic errors. 



2. 3. Controlled displacement by fast voltage switching 

The displacement operation |n) — > D(a)\n) = \a,n) is performed by switching the 
voltage on the neighboring segment to the trap segment, with its center located 
280 pm away from the initial trap position. A voltage of 1.0 V at this segment 
will give rise to an electric field of 600 V/m at the initial ion location, which is 
derived from electrostatic simulation of the trap [26J. This field exerts an acceleration 
of approximately 1.5-10 9 m/s 2 . Assuming an ideal, sudden switching process, the 



CONTENTS 



6 



potential minimum is shifted by a well defined distance along the trap axis for well 
defined time. The relation between the set voltage and the magnitude of the shift of 
the potential minimum is established via electrostatic simulation of the trap and has 
been experimentally confirmed [27] . For the displacement, the set voltage is changed 
instantaneously, however H-type filters with cutoff frequencies of approximately 300 kHz 
lead to smoothened voltage ramps at the trap segment. The actual displacement 
amplitude after the voltage kick is therefore lower than the value corresponding to 
the position space shift of the potential minimum, furthermore the displacement pulse 
is delayed and blurred. With the deviation of the potential minimum from its original 
location Xo(t), the resulting displacement is obtained from classical harmonic oscillator 
dynamics [20J: 

The fast switching of the voltages is done with a homebuilt fast multichannel 
arbitrary waveform generator. Its main building blocks are quad 16-bit serial digital 
to analog converters [|J which cover the output voltage range of -10 V to +10 V 
with a resolution of about 0.3 mV. The digital data is supplied in real time by a 
field programmable gate array (FPGA) [§J which in turn is programmed via a Gigabit 
Ethernet connection from a control computer. The maximum update rate is limited by 
the clock rate of the FPGA and is given by 2.5 MSamples/s, and a number of up to 
48 output channels can be controlled individually. A crucial feature is the possibility 
to vary the duration of each individual sample in steps of 20 ns beyond the minimum 
value of 400 ns, which allows for resolving the trap oscillation period [21] . 
Each output channel at voltage Vi is fed into a stage where a differential voltage Vd is 
added (subtracted) such that two channels with output voltages Vi ±Vd are obtained, 
which are fed to adjacent segments at the same axial location in the trap, which allow 
for the compensation of micromotion. The differential voltage Vd can be derived from a 
DAC, or from an external input. 

The measured output voltage noise at the trap frequency, in the steady state, is as low 
as 1.0-10~ 9 V/V Hz. This results in a spectral density of the electric field noise of about 
3.6- 10 -13 V 2 / Hz m 2 , which in turn yields a heating rate that lower than the actually 
observed heating rate of n w 0.3 ms _1 [24j by one order of magnitude; it therefore does 
not significantly contribute to the anomalous heating. Care was take to suppress digital 
noise, i.e. feedthrough of the FPGA clock and other digital signals onto the output 
channels is reduced by using high speed magnetic insulators for all digital signal lines 
and careful separation of analog and digital signals on the circuit boards bearing the 
DACs. 

We predict the resulting displacement magnitudes from simulations, which take into 
account the actual voltage waveforms applied to the shift segment. These waveforms are 
measured behind the filter. Averaging of about 10 5 identically generated waveforms was 

| DAC8814, Texas Instruments 

§ Avnet Virtex-5 FXT Evaluation Kit 
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necessary to make the results independent of measurement artifacts. These waveforms 
are used in conjunction with the simulated electrostatic potentials for the numerical 
solution of the classical equation of motion, 

rr dV A V) rr ^ 
mx = -e U A e U B {t) — f— , (5) 

where the forces are given by the spatial derivatives of the normalized electrostatic 
potentials v} 1V \x) along the trap axis, which result from a segment voltage of 1 V 
applied to electrode i. These forces are weighted by the respective voltages applied to 
the electrodes and added. The final energy Ef at times larger than the settling time 
is rescaled to harmonic oscillator units, Ef = Ef/hu, and the displacement is obtained 
from \a\ = yJ~Ef~. The resulting displacement amplitudes versus the kick amplitude Vk 
are shown in Fig. [3] 



3. Results 



Theory Measurement 




Figure 2. Resulting theoretical PPD a) and PPD reconstructed from the 
measurement result b) for the preparation of n = 1. The PPD are shown for each 
applied kick voltage. The interference minimum is clearly observable as a ditch bending 
toward higher phonon numbers for larger kick voltages. 



With the optically driven ladder climbing method described above, we prepare 
number states for n — 0, 1,2 before applying the kick. The displacement is determined 
by the magnitude of the voltage kick on the neighboring segment. The kick voltage 14 



CONTENTS 



8 




Figure 3. Extracted values for the displacement parameter a versus the kick voltage 
Vk for the datasets pertaining to n = (black) , n = 1 (blue) and n = 2 (red) . The solid 
lines are fits to polynomials of fourth order. The displacement parameters are obtained 
by fitting the theoretical predicition Eq. [7] to the reconstructed PPDs. The curves are 
clearly consistent, providing justification for the evaluation method. A slight deviation 
is observed for the n = case, which we attribute to a slow drift of the trap voltages. 
The green points show results from simulations, see text. 



is changed in steps of 200 mV, reaching a maximum displacement of a « 2.1 at Vk = 2 
V. In each case, the state reconstruction is performed with freely varying PPD. As a 
result, we obtain a normalized PPD for each pair of n and a. The results for n = 1 
are shown in Fig. [2j The measured PPD is compared to the theoretically expected one 

OS], 

pjW = |(A;|«,n)| 2 = e~\ a \ 2 \a\ 2{k - n)) n\k\ x 

We take into the imperfect preparation of the desired number state prior to the 
displacement operation. We assume that the action of the displacement kick is not 
dependent on the quantum state of the motional mode, which is justified by the 
spatial homogeneity of the electrical force. If we characterize the motional mode after 
preparation of the desired Fock state n by the density matrix pQ , we obtain the PPD 
which is actually measured from a convolution of the pure PPD from Eq. [6] with the the 
diagonal elements of the density matrix corresponding to the measurement with zero 
kick voltage, i.e. the PPD without displacement operation: 

^max 

Pk = Yj Tr (po n) \ m ) H) Pk ■ (7) 

m=0 

From the reconstructed density matrices without kick, we determine the respective 
Fock state preparation fidelities Tr (pj 0) |0) (0|) « 0.92, Tr (p^ |1) w 0.77, and 



" (_iy| a |2(-0 
t l\(n-l)\(k-l)\ 



(6) 



CONTENTS 



9 




1 2 30 1 2 30 1 2 

Displacement a 



Figure 4. Resulting PPDs for different number state preparations a) n = 0, b) n = 1, 
c) rt = 2. The respective PPDs are plotted versus the displacement amplitude a, which 
in turn is extracted by fitting the PPD to Eq. [6] The solid lines display the theoretically 
expected values from the same Eq. [6j The matching of the experimental values to 
the theoretical values becomes worse for larger displacement amplitudes because the 
precision of the maximum likelihood reconstruction is deteriorated by truncation errors. 
The insets for the n — 2 show bar charts of the PPDs at a « 1.5 and a ~ 1.5 to clearly 
show the interference minima. 



Tr(p}> ' |1) (1|J ss 0.72. To prove the consistency of this evaluation method, we show the 
extracted displacement magnitudes a versus the kick voltage for the three different 
data sets in Fig. |3j When these values are compared to displacement values inferred 
from simulations, we find an agreement on the single quantum level. 
In Fig. [2j we show the resulting PPDs for the case n = 1 plotted versus the kick 
voltage. The results clearly exhibit a single minimum in the PPD, which is moving to 
higher phonon numbers for larger kick voltages. The occurrence of two minima in the 
PPD is visible for n = 2 (see Fig. [4]), which has not been observed before in photonic 
realizations of DNS. The fact that n minima can occur in the PPD of a state \a,n) 
can be directly read off from Eq. [6j Within the sum a polynomial of order n in \a\ 2 
with coefficients of alternating signs is formed, which maximally exhibits n zeros upon 
squaring [9]. 

To provide a semiclassical explanation for the occurrence of the minima, we illustrate the 
idea of interference in phase space in Fig. [5] [2H]- The overlap integral (k\at, n) is given 
by the coherent sum of the two intersection areas A between the bands representing the 
number states \k) and |n). In the summation, one has to take into account the relative 
phase <fink{oi) — B/H which is given by the area B enclosed between the two classical 
trajectories: 

2 

\(k\a,n)\ 2 = — cos(t) nk a. (8) 

7171 

As the preparation number n increases, the values of <p n k{cx) sweep over larger ranges 
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Figure 5. Illustration how interference in phase space leads to the observed minima in 
the PPDs: We exemplify the overlap integral (3|a, 2). The Fock states are represented 
semiclassically as bands in phase space, spread around the classical trajectories 
corresponding to the energy eigenvalues (n + l/2)%v (dashed lines). The analysis 
number state basis is represented by bands centered around the origin, the DNS states 
\a,n) are displaced by the amplitude a. The overlap integral is determined not only by 
the diamond shaped red areas A, but also by the intersection area B of the trajectories 
(shaded in light red), see text. 



for increasing analysis quantum numbers k and increasing a, thus allowing for a larger 
number of interference minima. 

4. Outlook and conclusion 

We have demonstrated the creation of DNS \a,n) with number states of up to n = 2 
and a ~ 2.9, going significantly beyond the current state of the art of photonic 
experiments done on DNS so far. Emphasis is put on the fact that the displacement 
is created by electrical kicks, not by resonant excitation by rf pulses. This offers 
the flexibility to combine long-distance shuttling operations and controlled motional 
excitation on the single quantum accuracy level, which might open up the route to 
scalable quantum simulation schemes. The displacement operations were shown to be 
very accurate and reproducible. We emphasize that we were able to match the measured 
displacement amplitudes to simulation results on the single quantum accuracy level, 
where technological details such as the electrostatic properties of the microtrap and the 
voltage waveforms have to be determined with a sufficient precision. 
Future applications in the field of quantum state engineering include the joint usage of 
spin dependent optical forces and electrical forces, e.g. for optically assisted splitting 
of ion crystals or for phase referencing optical beat patterns. The latter might open up 
the possibility for novel measurement schemes where an optical interference pattern is 
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to be applied with a fixed phase relation to ion strings [29], which could be of use for 
quantum simulation experiments [30J. 
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